N=64;
j=sqrt(-1);
a=ones(N,1);
row=1;
for i1=1:N;
	a(i1) = sqrt(row+i1-1+.77777)+j*((row+i1-1+.22222)*(row+i1-1+.22222) / N - N/2);
end;
a(1) = N+3+j*imag(a(1));
a(2) = 1-N+j*imag(a(2));
A=fft(a);
a1=ifft(A);
if N<=32
	A
	a1
end
max(abs(A-atrans))./(max(abs(A))/N)
max(abs(a-aitrans))./(max(abs(a))/N)
%max(abs(a1-aitrans))./(max(abs(a1))/N)
